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Résumé 

Dans ce travail, nous étudions un système non-linéaire structuré en âge et maturité 
décrivant la production des cellules sanguines dans la moëlle osseuse. Après intégration 
du modèle, nous obtenons une équation aux dérivées partielles du premier ordre faisant 
apparaître un retard distribué en temps et une dépendence non-locale dans la variable 
maturité. Nous prouvons que l'unicité des solutions ne dépend que des cellules de petites 
maturités (cellules souches) et nous donnons un résultat d'invariance. 

1 Introduction 

Nous nous intéressons dans ce travail à un modèle mathématique non-linéaire décrivant 
une population cellulaire : la production des cellules sanguines dans la moëlle osseuse, à 
partir des cellules souches. Ce modèle est décrit par un système de deux équations aux 
dérivées partielles structurées en âge et maturité. 

Les modèles de populations biologiques structurés en âge et maturité sont apparus vers la 
fin des années 60. Les premiers à les avoir étudiés sont Keyfitz [Ij (1968), PoUard [S] (1973), 
Henry [3] (1976) et Mackey [ïïj (1978). Le modèle de Mackey, dans lequel la population 
cellulaire est supposée uniquement proliférer, a été étudié plus amplement à partir de 1990, 
notamment numériquement par Mackey et Rey dans [TD] (1992) et [TT] (1993). En 1994, 
Mackey et Rudnicki [7] ont complété ce modèle en tenant compte d'une phase de repos 
en plus de la phase de prolifération. Ils obtiennent les équations sur lesquelles nous allons 
travailler et qui sont présentées en l2.1l et l2.2l Dans la phase de prolifération, dont la durée est 
limitée, les cellules synthétisent de l'ADN puis se divisent, donnant naissance à deux cellules 
filles. Chacune des cellules filles entre ensuite immédiatement dans la phase de repos. Elle 
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peut y demeurer toute sa vie et y mourir de vieillesse ou bien être introduite dans la phase 
de prolifération à un certain moment. Nous étudirons ici le cas non-linéaire où le taux de 
réintroduction dépend de la population totale de cellules au repos. 

Mackey et Rudnicki montrent que le système d'équations aux dérivées partielles qu'ils ont 
obtenu devient, après intégration, un système d'équations aux dérivées partielles à retard. 
Dyson, Villella-Bressan et Webb en 1996 [2] et Adimy et Pujo-Menjouet en 2001 [T] ont 
montré que la population de cellules de petites maturités influence le comportement de la 
population totale, dans le cas où les cellules sont supposées se diviser toujours au bout du 
même temps. 

Nous nous proposons d'étudier le modèle proposé par Mackey et Rudnicki en 1994 [7] 
en tenant compte d'une remarque faite par Mackey et Rey dans |11| concernant la division 
cellulaire : nous supposerons que l'âge de la cytocinèse (le point de division cellulaire) est 
distribué selon une densité de support [t, r], avec < r < r < +oo. Dans [TT], Mackey et 
Rey n'ont étudié que numériquement ce modèle et sans tenir compte de la phase de repos. 
Dans les travaux de Mackey et Rudnicki [7], Dyson, Villella-Bressan et Webb [2] et Adimy 
et Pujo-Menjouet [T], les cellules proliférantes se divisaient toujours au bout du même temps 
fixe. 

Enfin, rappelons que la production de cellules sanguines peut avoir deux comportements 
différents. Le premier correspond à une production normale de cellules et a lieu lorsque la 
population de cellules souches (c'est-à-dire les cellules de petites maturités) est strictement 
positive. Le deuxième mène à une production anormale de cellules et est le résultat de la 
destruction ou de l'absence de cellules souches dans le milieu originel : ce cas correspond à 
l'anémie aplasique, une maladie s'attaquant aux cellules souches. 

En supposant que la phase de prolifération est suffisament longue, nous allons montrer que 
l'unicité des solutions dépend uniquement des cellules souches et nous en déduirons un résultat 
d'invariance soulignant l'influence des cellules de petites maturités sur le comportement de 
la population totale. 

Notre travail est organisé de la façon suivante : dans la section [2l nous présentons le 
modèle biologique sur lequel est basé notre travail ; dans la section [3l nous obtenons, par 
intégration, le système d'équations aux dérivées partielles décrivant l'évolution de la popula- 
tion; dans la section 31 nous donnons une formulation intégrée du problème; dans la section 
O nous établissons un résultat d'unicité des solutions ne dépendant que des cellules de petites 
maturités et dans la section [5] nous donnons un résultat d'invariance. 

2 Présentation du modèle 

La vie d'une cellule débute avec la division d'une autre cellule, dite cellule mère, et se 
termine, généralement, lors de sa propre division; la cellule donne alors naissance à deux 
nouvelles cellules appelées cellules filles. Durant sa vie, chaque cellule suit le même chemin 
que celui pris par sa mère : il s'agit du cycle cellulaire. Le cycle cellulaire tel que nous 
l'entendons ici comprend deux phases : une phase de prolifération et une phase de repos. 

La phase de prolifération est la phase active du développement cellulaire : elle est dédiée 
à la synthèse de l'ADN et à la division cellulaire. Elle s'achève lors de la cytocynèse : la 
division de la cellule proprement dite. Une cellule ne peut rester indéfiniment dans la phase 
de prolifération : si elle ne se divise pas, elle meurt avec un taux 7. 
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Juste après leur naissance, les cellules filles entrent immédiatement dans la phase de repos, 
ou phase Gq (de l'anglais Gap), qui est un stade quiescent du développement cellulaire : les 
cellules au repos sont inactives. Contrairement à ce qui se passe dans la phase de prolifération, 
les cellules peuvent rester toute leur vie dans la phase de repos. La plupart seront introduites 
à un certain moment dans la phase de prolifération et compléteront ainsi le cycle cellulaire, 
mais les cellules au repos peuvent aussi mourir, avec un taux ô. 

Une cellule sera caractérisée par son âge et sa maturité. La maturité désigne un aspect 
morphologique du développement de la cellule, au même titre que la taille ou l'âge ; il s'agit 
d'une variable continue : dans l'évaluation de la maturité, on tient compte, par exemple, du 
niveau de synthèse d'ADN, d'ARN, de la présence de mitochondrics, etc. La définition de la 
maturité est liée à la nature de la cellule. Nous supposerons que la maturité m d'une cellule 
varie de m = à m = 1. 

Toutes les cellules proviennent des cellules embryonnaires issues des premiers joins du 
développement, appelées cellules souches {stem cells). Il s'agit de cellules indifférenciées et 
immatures (c'est-à-dire de petites maturités) et capables, d'une part, de se multiplier (et 
donc d'augmenter la quantité de cellules souches), et d'autre part, de donner naissance à des 
cellules différenciées (cellules du foie, cellules musculaires, cellules sanguines, . . .). 

Il existe quatre types de cellules souches dont les cellules souches hématopoïétiques de la 
moëlle osseuse. Ce sont celles-ci qui nous intéressent car elles sont à l'origine de toutes les 
cellules sanguines : globules rouges, blancs, plaquettes, etc. 

Mis à part leur importance dans le développement d'une population cellulaire, il faut 
noter que les cellules souches sont aussi à l'origine de certaines maladies, graves, dont l'anémie 
aplasique. Il s'agit d'une maladie rare et sérieuse caractérisée par l'absence de formation ou 
de division des cellules souches hématopoïétiques, entraînant une chute de la production des 
cellules sanguines. 

Plusieurs facteurs expliquent l'apparition de cette maladie. Dans 50 à 65 % des cas la 
nature de la maladie est idiopathique ou primitive, c'est-à-dire qu'elle ne comporte pas de 
cause connue. Dans les autres cas, elle est héréditaire ou duc notamment à des expositions aux 
radiations nucléaires, à des toxines de l'environnement (insecticides, vapeurs de benzène, . . . ) 
ou bien au traitement du cancer par radiothérapie, chimiothérapie ou d'autres médicaments. 

2.1 La phase de repos 

Nous désignons par n{t, m, a) la densité de cellules dans la phase de repos à l'instant t, 
de maturité m et d'âge a. Nous notons 




'OO 



n{t, m, a)da, 



le nombre total de cellules de maturité m dans la phase de repos à l'instant t. 
La densité de cellules n{t, m, a) vérifie l'équation de conservation : 



dn dn d{Vn) 



= -{S + (3)n, 



(1) 



dt da dm 



oii S est le taux de mortalité, /? le taux de réintroduction des cellules dans la phase de 
prolifération et V : [0, 1] — > [0, +00) est la vitesse de maturation des cellules. 
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5 = 5{m) et f3 = f3{m, N{t,m)) sont supposées positives et continues, et /3 est supposée 
décroissante par rapport à N{t,m). V est de classe sur [0,1] et vérifie V{0) = et 
y (m) > pour tout m G (0, 1]. De plus, nous supposerons que 



/ — -— = +00, pour tout m g (0, 1 



Cela signifie qu'une cellule de maturité nulle ne peut jamais atteindre un état où sa maturité 
serait non nulle. 

Exemple 2.1 Si V{s) as^ , pour tout s G [0, 1], avec a > et p > 1, alors Jq™ est 
une intégrale de Riemann divergente pour tout m G (0, 1]. 

n satisfait de plus la condition initiale : 

n(0, m, a) — T(m, a), pour tout (m, a) G [0, 1] x [0, +00), (2) 

où T est une fonction suffisamcnt régulière vérifiant lima^oo T(m, a) = 0. 

2.2 La phase de prolifération 

Dans la phase de prolifération, l'âge d'une cellule varie de à t où t est la durée maximale 
de cette phase : parvenue à l'âge r une cellule doit s'être divisée sinon elle meurt. Cela entraîne 
qu'une cellule ne peut rester indéfiniment dans la phase de prolifération. 

Soit p{t, m, a) la densité de cellules de maturité m et d'âge a dans la phase de prolifération 
à l'instant t. Nous notons 



P{t,m)= f 



p(t, 771, a)da, 



le nombre total de cellules de maturité m dans la phase de prolifération à l'instant t. 

Soit 7 une fonction de la maturité m, positive et continue, représentant le taux de mortalité 
des cellules dans la phase de prolifération. 

p{t, m, a) vérifie l'équation de conservation : 

dp dp d{Vp) 
ot oa dm 

où V représente la vitesse de maturation des cellules dans cette phase, définie en l2.1l Nous 
supposons que la vitesse de maturation des cellules est la même dans les deux phases. 
p vérifie en plus la condition initiale : 

p{0,m,a) — T{m,a), pour tout (m, a) G [0, 1] x [0, t], (4) 

où r est une fonction suffisament régulière, définie sur [0, 1] x [0,t]. 

Pour être bien posé, le système d'équations (IT|)-(l3|) doit être muni de conditions aux bords, 
traduisant le comportement des cellules lors des changements de phases. 
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2.3 Conditions aux bords 

En accord avec les remarques de Mackey et Rey ([H]), nous supposerons que le point de 
cytocinèse des cellules proliférantes est distribué sur un intervalle [r, t], oii < r < r < +00, 
selon une densité a 1— > k(rn,a). r représente la durée maximale de la phase de prolifération 
tandis que t est choisi strictement positif car les cellules ne se divisent immédiatement après 
être entrées en phase de prolifération que dans le cas de certaines maladies qui ne nous 
intéressent pas ici. 

Nous notons g{m) la maturité d'une cellule fille à sa naissance, lorsque la maturité de 
la cellule mère valait m. La fonction g : [0, 1] [0, 1] est supposée continue, strictement 
croissante sur [0, 1], de classe sur [0, 1) et telle que g{m) < m pour m & {0,1). g est alors 
inversible et nous prolongeons g~^ en posant g^^{m) — 1 si m > g{l). 

L'application k est telle que m 1-^ k{m, a) est la proportion de cellules mères d'âge a dans 
la phase de prolifération pouvant donner après division une cellule fille de maturité m. k est 
supposée positive et continue sur [0, 1] x [r, r], et telle que k{m, a) = si m > g{l)- Ceci par 
définition de la fonction g : la maturité d'une cellule mère ne pouvant être supérieure à 1, 
celle d'une cellule fille ne peut excéder la valeur 5(1). Nous considérerons donc, dans toute 
la suite, que la maturité d'une cellule fille est comprise entre et 5(1). 

La première condition au bord est alors donnée par 

n{t,m,0) = 2 / k{m,a)p{t,g~^{m),a)da. (5) 



Elle signifie qu'une cellule de maturité m peut provenir d'un certain nombre de cellules 
mères, d'âges différents, contrairement aux modèles précédemment étudiés par Mackey et 
Rudnicki [7], Dyson, Villella-Bressan et Webb [5] et Adimy et Pujo-Menjouet [T], où la division 
survenait toujours au bout d'un temps fixe (la fin de la phase de prolifération). 

La seconde condition aux limites est 



00 



p{t,m,0)^ j3{m,N{t,m))n{t,m,a)da = P{m,N{t,m))N{t,m). (6) 





Elle représente le flux de la population de cellules, de la phase de repos vers la phase de 
prolifération. 

Nous considérons le flot tTj : [0, f ] [0,1], défini pour s < 0, solution de l'équation 
différentielle ordinaire : 

[ u{0) = m, 

qui représente l'évolution de la maturité des cellules pour atteindre une maturité m à l'instant 
à partir d'un temps s < 0. tTs vérifie noijn) = m, 7rs(0) = et Trs{m) g (0, 1] pour s < et 
me (0,1]. 

En fait, nous pouvons remarquer (cf. [9]) que tts peut s'écrire explicitement 
TTsim) = h~^{h{m)e^), pour m G [0, 1] et s < 0, 



011 h : [0, 1] [0, 1] est défini par 




pour m S (0, 1], 
pour m — 0. 
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Notons que h est une fonction continue et strictement croissante, donc inversible. 
Exemple 2.2 Si V{m) — am, pour tout ni G [0, 1], avec a > 0, alors 

h{m)^m^l°', Vme[0,l]. 

Donc, pour tout m G [0, 1], 

h~^{m) = m", 

et 

7rs(m) — me"" , pour s < 0. 

Nous obtenons les fonctions utilisées par Rey et Mackey dans flO}/ et Jll}/ et Dyson, Villella- 
Bressan et Webb dans 0/. 

3 Equations du modèle 

Avant de déterminer les équations du modèle sur lesquelles nous allons travailler nous 
établissons un lemme qui nous sera utile. 

Soit A l'application définie de [0, +oo) x [0,g(l)] dans M+ par 

A(s, m) = 7r_,(.g-i(TO)) ^ h-\h{g-\m))e-''). 

Lemme 3.1 A satisfaits les assertions suivantes : 

(i) l'application s G [0, +oo) ^ A(s,m) est décroissante pour chaque m G [0,17(1)], 

(ii) l'application m G [0,(7(1)] A(s,m) est croissante pour chaque s > 0, donc, pour 
tout s > ei tout m G [0,5(1)], 

= A(s,0) < A(s,m) < A(s,<7(l)) = h~\e-"), 

(m) pour tout m G (0,(7(1)], 



A(s,m) < m ssi s > 
(iv) pour tous a,s>0, et tout m G [0,.g(l)], 



Vie)' 



7r_CT(A(s, m)) — A(s + a, m). 

Preuve : Les points (i) et (ii) sont évidents par définition, les fonctions h et g étant crois- 
santes sur [0, 1]. 

, . / , /-s^'M de , h{g-\m)) 

Pour montrer le pomt (m), remarquons que / — r-r- = In — -- — 7 — pour tout 

7™ V{e) h{m) 

m > 0. Alors nous avons, pour tout m G (0,(7(1)] et tout s > 0, 

A(s,to) < m <^ h{g^^{m))e^" < h[m) 

h(jn) 



^ -s < In 



<^ s > 



h{g-\m)) 
_M_ 

Wr 
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D'où le résultat. 

Montrons (iv). Par définition nous avons, pour tous cr, s > et tout m G [0, (/(!)], 

TT-a{A{s,m)) = h-\h{h'\hig-\m))e-''))e-'') 
= h-^{h{g-^{m))e-^'+^^) 
= A{s + a,m). 

Le lemme est ainsi démontré. □ 

Nous utiliserons, dans toute la suite, les notations suivantes : 

^{m,t) = ea;p| - J 7(7r_s(m)) + F'(7r_s(m))ds|, 

pour tout m G [0, 1] et tout t > 0, et 

C(m,a) = k{m,a)^{g~'^{m),a), 

pour tout m G [0, 1] et tout a > 0. 
Nous établissons le résultat suivant. 

Proposition 3.1 Le nombre total, N{t,m), de cellules de maturité m dans la phase de repos, 

vérifie les équations suivantes : 

(i) SiO<t<Tetme [0,5(1)], 

-N{t, m) + —{V{m)N{t, m)) = -{5{m) + /3(m, N{t, m)))N{t, m) 

+ 2^{g-'^{m),t) J k{m,a)T{A{t,m),a-t)da. (7) 

(ii) Siz<t<TetmG [0,5(1)], 

-N{t, m) + —{Vim)N{t, m)) = -(5(m) + /3(m, N{t, m)))N{t, m) 

+2 J C{m,a)P{A{a,m),N{t - a,A{a,m)))N{t - a,A{a,m,))da 

+ 2^{g-^{m), t) k{m, a)T{A{t, m), a - t)da. 

(m) Sit>T etm.e [0,5(1)], 
d d 

Q-^N{t,m) + —{V{m)N{t,m)) = -{S{m) + P{m,N{t,m)))N{t,m) 



(8) 



+ 2 / C{m,a)l3{A{a,m),N{t - a,A{a,m)))N{t - a,A{a,m))da. (9) 
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(iv) Sit>0 etme [g{l), 1], 

— iVU, m) + —(V(m)N(t, m)) = -(S(m) + Mm, NU, m)))N{t, m). (10) 
at dm 

De même, la population totale, P{t,m), de cellules de maturité m dans la phase de pro- 
lifération vérifie les équations suivantes : 

(v) SiO<t<T etme [0, l], 

d d 

—P(t,m) + —(V(m)P(t,m)) = -'y(m)P{t,m) + (3{m, N(t,m))N(t,m) 
ot dm 

- T{Ti^t{Tn),T-t)^{m,t). (11) 

(vi) Si t >T et m [0, 1], 

d c) 

^P(t, m) + —(V(m)P{t, m)) = -j(m)P(t, m) + (3{m, N(t, m))N(t, m) 
ot dm 

- ^im,T)f3{TT_-{m),N{t - t, TT.-{m)))N{t - t, n_-{m)). (12) 

Preuve : Elle est basée sur la méthode des caractéristiques. Nous intégrons les équations JT]) 
et ([3]) par rapport à l'âge : 

d d 

—P(t,m) + ^—(V(m)P(t,m)) = -j(m)P(t,m) + p(t,m,0) 
ot dm 

-p{t,m,T), 

d d 

— N(t,m) + —(V(m)N(t,m)) ^ -(ô(m) + l3(m, N(t,m)))N(t,m) 
ot dm 

+n{t, m, 0), 

en ayant supposé que lima^oo 'n{t, m, a) = 0, ce qui signifie simplement, d'un point de vue 
biologique, que toutes les cellules sont condamnées à mourir. 
En utilisant la méthode des caractéristiques nous obtenons 

p(t,™,a) = (^J^,,r_.(m), a-t)Ç(m,t) si ^ t < a, 
[p(t — a, T:-a(m),0)^(m, a), SI t > a. 

Nous pouvons alors écrire, en utilisant 

p{t,m,T) = p{0,TT^t{m),T - t)S,{m,t), 
= r(Tr^t{m),T - t)^{m,t), 

si < t < T, et 

p{t,m,T) = p{t -T,Tr^r{m),0)^{m,T), 

= P{T:Mm),N{t - T, T:^-{m)))N{t - T, 7r_-(m))e(m, t), 

si i > r. 
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En utilisant la relation ([S]), nous obtenons bien que P{t,m) vérifie les relations (|lip et 

M- 

De même, d'après ^ et en utilisant (|13p . nous pouvons écrire : 



n(t,TO,0) ~ ^ J ^(™jfl)p(^:5 ^{m),a)da 

k{m,a)p{0,A{t,m),a - t)^{g^'^{m),t)da, 
= 2^{g-\m),t) [ k{m,a)T{A{t,m),a-t)da, 



si < t < r, et 

et 



n(t,m,0) ~ ^ J ^(™;o)p(i,(7 ^{m),a)da 

+2^ k{m,a)p{t,g~^{m),a)da, 
~ J k{m,a)p{t — a,A{a,m),0)^{g^^{in),a)da 

+2^ k{m,a)p{0,A{t,m),a- t)^{g^'^ {■m),t)da, 
= '2 J k{Tn,a)f3{A{a,m),N{t - a,A{a,m)))N{t ~ a,A{a,m)) X 

C{g~'^im),a)da + 2^{g^^{m),t) k{m, a)T{A{t, m), a - t)da, 

si T < t < T, et 

n(t,m,0) — 2 k{m,a)p{t a,A{a,m),0)^{g^^{m),a)da, 



= 2y k{m,a)P{A{a,m),N{t ~ a,A{a,m))) X 
N{t - a, A(a, m))£,{g~^ (m), a)da. 

S\t>T. 

N{t,m) vérifie donc bien les équations ^ et 

En remarquant que, pour tout m G [^(l), 1], Cl"';'^) = et k{m,a) — 0, nous obtenons 
l'équation (fTÏÏ| .a 



Nous remarquons tout d'abord que les solutions des équations ([7]), ([5]) et pT|) deviennent 
des conditions initiales pour les équations ([9]) et respectivement. De plus, les solutions de 
l'équation ([9]) ne dépendent pas des solutions de l'équation (fT2)) . Enfin, si nous connaissons le 
comportement des solutions de ^ nous pouvons en déduire les solutions de l'équation (fT2|) 
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par la méthode des caractéristiques. Nous concentrerons donc notre étude sur les équations 
décrivant le comportement des cellules dans la phase de repos. 

Nous allons déterminer une formulation faible du problème ©-([SI)-®. 

4 Formulation faible 

Soit C*' [0,(7(1)] l'espace des fonctions continues sur [0,17(1)], muni de la norme du maxi- 
mum, et soit A : X'(A) C C°[0,g(l)] C°[0,5(l)] l'opérateur défini sur l'ensemble : 

V{A) = (u G C"[0,g(l)];u G C^{Q,g{l)], lim V{m)u'{vi) = o|, 

par 

_ r -{5{m) + V'{m))u{m) - V{m)u'{m), si m e (0,5(1)], 
^^"\-((5(0) + r(0))w(0), sim = 0. 

Proposition 4.1 A est le générateur infinitésimal du C^-semi-groupe {S(t))t>o défini, pour 
ip G 2?(A), m G [0,5(1)] ett>0, par : 

{S{t)iP){m) = iP{h-\h{m)e-*))exp!^ ~ ((5 + ^')(^"^(^Me"'))rfs|- 



Remarque 4.1 Comme h ^{h(m)e^) — 7rs{m) pour tout m G [0,1] et tout s < 0, alors, 
\fip G V{A), Vm G [0,5(1)] etyt>0 : 

(5(0V')(™) - ^(^^t(m))if (t, m), (14) 

où ^ 

K{t,m)=exp\ - / ô{Tr_sim)) + V'{TT-s{m))ds\, 



pour tout m G [0,5(1)] et tout t >0. 

Preuve de la Proposition : La démonstration de cette proposition nécessite l'utilisation 
de trois Icmnics. 

Soit B : V{B) C C°[0,5(l)] ^ C°[0,5(l)] l'opérateur défini par : 

V{B) = {u G C"[0,5(l)] ;u G C^{Q,g{l)], lim V{m)u{m) = o|, 

Bu{m) = I ^^M«'M' si ^ e (O'ff(l)]' (15) 
^ ' [ 0, SI m = 0. 

Lemme 4.1 L'opérateur B défini par hl5\) sur D{B) est le générateur infinitésimal d'un 
semi-groupe de contractions sur C'^[0, g(l)]. 
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Preuve : Nous allons utiliser le théorème de Hille-Yosida. B est un opérateur linéaire fermé. 
De plus, en remarquant que C^[0,(7(l)] C 2?(iî) C C°[0,(ji(l)] et comme C^[0,g(l)] est dense 
dans C°[0, (/(l)], nous obtenons bien que Î?(-B) est dense dans C°[0,(7(l)]. 
Soient w G C''[0,(7(l)] et A > donnés. Montrons qu'il existe une fonction u G 2?(S) et une 
seule telle que 

{I ~ \B)u = w, 

011 / est l'opérateur identité, et que 
011 la norme |.|oo est définie par 

\u\oo — sup \u{m)\. 

0<m<g(l) 

Si u existe alors u est solution de l'équation différentielle 



u(m) + XV{m)u'{m) — w{m), pour < m < g{l), 
u{0) = w{0). 



(16) 



Considérons la fonction h\ définie par 

ds 



\ exp( ), pour m € (0,g(l)], 

(to) = S ^7g(l)AF(s)y 

0, pour m = 0. 



Les solutions de l'équation p6|) sont de la forme 

h\{m) \ JgÇT^^ XV{s) J 
pour m G (0,g(l)] et fc G M. h\ est continue sur [0,g(l)], de classe sur (0,g(l)] et 

, feA(™) 

pour m G (0,g(l)]. Donc 

Mf)d. = /.,(m)-/i,(g(l)). 

9(1) ^^(s) 

Ainsi, lim„i^o xp^^^s = -'iA(g(l)) = -1- Par conséquent, G 2.1(0,5(1)) et, 

comme /ia(0) = 0, la fonction u donnée par p7p est bornée si et seulement si 

Ân;)'^^^^)''" 

L'unique solution possible de l'équation (jl6p est donc 



h\(m) \Jq XV{s) 
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pour m e (0,.g(l)], et u{0) = w{0). 

De plus, si M G T^{B), alors \u{m)\ < |w|oo pour tout m G (0,^(1)] et, comme u{0) = w(0), 
\u\oo < \w\oo- 

En appliquant le théorème de Hille-Yosida, nous achevons la démonstration du lemme. □ 



Lemme 4.2 Si ip € T>{B), alors 

lim — — i — - = -Vlmm (m). 



Preuve : Nous effectuons le changement de variables suivant : 

y = h-\h{m)e-'). 

Alors, par définition de h, 



t = In 

Ainsi 



ds 



h(rn) \ 



lim '^^^ ^{h{m)e — ip{m) ^. il){y) — 'ip{m) y — m 



t^o+ t v^m y — m — 



y v(s) 
= -y(m)V'(m). 



Ce qui conclut la démonstration. □ 



Par définition du générateur infinitésimal d'un semi-groupe, nous déduisons de ce lemme 
que le semi-groupe engendré par B, que nous notons {T{t))t>o, est donné, pour tout tp € 'D{B) 
et tout t > 0, par 

(T(t)V')(m) = ^(/î-i(/i(m)e-*)), Vm G [0,5(1)]. 

Lemme 4.3 Si {T{t))t>o est le semi-groupe engendré par B alors le semi-groupe engendré 
par B -\- al, avec a e ^^[0, g{l)], est défini par 

{S{t)ij){m) = exp(^ J {T{s)a){fn)ds^ {T{t)i^){m), 
pour tout ip G 'D{B), tout m G [0,5(1)] et tout t > 0. 

Preuve : Il suffit de montrer que pour tj) G T>{B), 



(S'(t)V')(m) - V(m) ,„ X 

hm ^ ^ '^^^ ' = (B-\- al)é(m). 

t^o+ t ^ /rv / 



12 



M. Adimy, F. Crauste 



Un modèle de prolifération cellulaire 



D'après la définition du semi-groupe {S{t)) nous avons : 

(5(t)^)(m)-V.(m) ^ lï^ /\T(.)a)Md.)(TWV.)M-V.M 



V "'0 

exp[ / {T{s)a){m)ds^ ' ^ ^^^'^ ' ' 







V t 



exp[ j {T{s)a){m)ds ]i/j{m) — ip{m) 



En passant à la limite sur t nous obtenons le résultat : le premier terme de droite tend vers 
B'ip{m) et le deuxième vers {al)^{ni). □ 

En utilisant les lemmes W?]\ et avec a = ô + V' , nous obtenons le résultat énoncé dans 
la proposition. □. 

Proposition 4.2 Soit lç g C"([0,t] x [0,(7(1)]). La formulation intégrée du problème 
(Ml-ÏB est : 

N{t,m) = ip{T,TT_(^t_y){m))K{t -T,m) + G{N){t,m) 

-J{N){t,m) sit>T, (18) 

N{t, m) = ip{t, m), si < t <T, 

pour tout m G [0,g(l)], avec 
GiN){t,m) = 2j_ ^ C(7r_(t-.)(™),a)/3(A(a,7r_(t_,)(m)),7V(s-a,A(a,^_(t„,)(m))) x 

N{s — a, A(a, Tr^çt-s) {m)))dci^ K{t — s, m)ds, 

et 

J{N){t,m) = K(t ~ s,m)/3(7r_(t_s)(™),^^(s,7r_(4_s)(m)))iV(s,7r„(4_s)(»Ti))ds, 
pour tous i > T ei m G [0, 5(1)]- 

Preuve : L'équation ([9]) se met sous la forme : 

d d 
—N{t,m) = -(ô(m) + l3{m,N(t,m)) + V (m))N{t,m) - V(m]—N(t,m] 
ot dm 



+ 2j C{m,a)P{A{a,m),N{t ~ a,A{a,m)))N{t - a,A{a,m))da. (19) 

Soit u l'application définie pour tout t > par : 

u{t){m) ^ N{t,m), Vtog[0,1]. (20) 
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Comme 

N{t - a, A(a, m)) = u{t - a)(A(a, m)) := u^ait){A{a, m)), 
il vient, d'après ([TÏÏ]), pour tout m G [0, 1] : 

u'{t){m) = -((5(m) + V'{m))u(t){m) - V{m){u{t)y (m) - (3{m, u(t){m))u{t){m) 

+2 / C{m,a)P{A{a,m),u{t - a)(A{a,m)))u^a(t){A{a,m))da, 



c'est-à-dire 

u'{t){m) = Au{t){m) - (3{m,u{t){m))u{t){m) 

+2 [ ({m,a)l3{A{a,m),u{t - a){A{a,m)))u-a{t){A{a,m))da. 



u{t) est donc donné, pour tout m G [0,(7(1)], par la formule : 

u{t){m) = S{t -T)ip{T){m) - S'(é - s)/3(7r_(4_s)(TO),u(s)(7r_(t_5)(m))) x 

u(s)(7r_(t„,)(TO))ds + 2 S{t-s)(^J^ C{m,a) x 
/3(A(a, m),u{t — a)(A(a, m)))it_a(s)(A(a, m))da ) ds, 



si t > r, et u{t){m) = tf{t){m) := ip{t,m) sit <t, où (p{t) G C''([0, g(l)]) pour tout t > 0. 
D'après p4)) . si ^ > r. 



u{t){m) = tp{T,n_(^t--){m))K{t - T,m) - (3{'K_^t-s){m),u{s){TT_(^t_^-^{m))) x 
u{s){TT_^t-s){m))K{t - s,m)ds + 2 (^j C{T^-{t-s){m),a) x 
/3(A(a, m), «(i — a)(A(a, m)))M_a(s)(A(a, 7r_(t_s) (m)))(ia^ K{t — s, m)ds. 
En utilisant (|2Ô1) nous obtenons alors (fTHl). □ 



5 Existence et unicité 

Nous allons à présent établir l'existence et l'unicité des solutions du problème ([T5)) . en 
utilisant une méthode itérative et la méthode des pas, avant de montrer que l'unicité des 
solutions ne dépend en fait que des cellules de petites maturités. 

Nous faisons l'hypothèse suivante : 

(Hi) Il existe / G R+ tel que Vto g [0, 5(1)], Vxi, a;2 G M+, 

\xi[3{m,xi) - X2P{m,X2)\ < l\xi ~ X2I. 
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Théorème 5.1 Sous l'hypothèse (Hi), pour toute donnée initiale ip G C''([0,t] x [0,.g(l)]), 
le problème fi<$|) possède une unique solution qui est continue. 



Preuve : Soit (ç G C"([0,t] x [0,(7(1)]) donnée. Si N est solution du problème (fT8|) alors 

N{t,m) = ip{T,TT_^t-r)im))K{t ^T,m) + G{N){t,m) - J{N){t,m), (21) 

pour tout t>Tei tout m G [0,.g(l)], et N{t,m) = (p{t,m) pour t G [0,t] et m G [0,g(l)]. 

Soit t G [0,r + t]. Nous définissons la suite (7V„)„gN de la façon suivante : 
si i G [0,r] alors 7V„(t, m) = (p{t,m), Vm G [0,.g(l)], et si i G [t,t + t] alors, Vm G [0,^(1)], 

No{t,m) = ip{T,Tr_i^t-T){nT'))K(t -T,m) + G{(f){t,m), 
N„{t,m) = Noit,m) - J(iY„_i)(i, m). 

Nous avons alors : 

\Ni{t,m) - Noit,m)\ = \J{No)it,m)\. 

Les fonctions S et V étant continues sur [0, .g(l)], l'application K est continue sur [0, +00) x 
[0,(?(1)] et donc bornée sur [r, r + r] x [0,5(1)]. Notons 



Alors, 



a = max l^(^i 1^)1 

(t,m)6[r,r+T]x[0,s(l)] 



et 

\Ni{t,m) ~ Na{t,m)\ < a 1 \/3{TT_^t-s){m), No{s,7T_çt-s){m)))\ x 



|iVo(s,7r_(t_5)(m))|ds, 
< al \No{s,n_(^t-s)im))\ds. 

(f, K et G{(p) étant continues, iVo est continue sur [t, r + r] x [0,.g(l)] : 

3M > 0, Vrn G [0,.g(l)], Vt G [t,t + t], \No{t,m)\ < M. 

Alors, 

\Ni{t,m) ~ No(t,m)\ <Mal(t-T). 
Montrons par récurrence le résultat suivant : 

(f _ 

(H„) ViG [r,r + T],VmG [0,g(l)],|7V„(f,m)-A^„_i(t,m)| <M(aO"^ r^- 

71 ! 

La relation est vraie pour n = 0, nous venons de le vérifier. Supposons (H„) satisfaite et 

montrons que (Hn+i) est réalisée. 

Soient t G [t,t + t] et m G [0,51(1)]. Alors, 

|iV„+i {t, m) - iV„ {t, m) 1 = 1 J(A^„) (i, m) - J(A^„_i) (t, m) | , 



15 



M. Adimy, F. Crauste 



Un modèle de prolifération cellulaire 



donc, en utilisant l'hypothèse (Hi), 

\Nn+i{t,m) - Nn{t,m)\ < ot^ j 7i'_(t_s)(m)) 

-7V„_i(s,7r_(4_s)(TO))|ds. 
En utilisant l'hypothèse de récurrence, nous obtenons 



|7V„+i(i,m)-iV„(t,m)| < M{alY+^ / -^ds, 

(f _ rV'+l 
< M(al)-+' ^) l> . 

(n + Ij! 

Nous obtenons bien le résultat escompté. 

Par conséquent, N := lim„^oo iV„ existe uniformément sur [t, r + t] et est continue. 
Montrons que N est solution de (HH]) pour t e [t, r + r] et m G [0,g(l)], c'est-à-dire 

N{t, m) = No{t, m) - J{N){t, m), 

car G{N){t,m) = G{(p){t,m), pour tout t G [t,t + t] et tout m G [0,(7(1)]. En effet, si 
a G [t, t] et s g [t, t] alors s — a G [0, t — r] C [0, r] et donc 

N{s — a, A(a, m)) = ip{s — a, A(a, m)), 

pour tout m G [0,g(l)]. Posons 

v{t,m) = |7V(t,m) - No{t,m) - J(A^)(t,m)|, 

pour tout t G \t,t + t] et tout m G [0,(7(1)] et montrons que v = 0. 
En remarquant que 

00 

N{t,m) - N„{t,m) = ^(iVp+i(t,m) - Np{t,m)), 

p—n 
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alors 



v{t,m) < \N{t,m) - N„+i{t,m)\ + \J{Nn){t,m) - J{N){t,m)\, 
< \N{t,m) - Nn+i{t,m)\ 



+al |A^„(s,7r_(i_s)(m)) ~ iV(s,7r„(t_5)(m))|ds, 

OO 

< J2 \Np+i{t,m)-Np{t,m)\ 



p— n+l 

et OO 



+al / ^ |iVp+i(s,7r_(t_s)(m)) - iVp(s,7r_(t_s)(m))|ds, 

p—n 



-n+l ^ ' p—n 



< 2M j2 (^^y^ 



ip 


+ 1)! 


[t- 




ip 


+ 1)! 


lit- 





-2 



-n+l p—n 



En faisant tendre n vers +oo, nous obtenons : 

v(t,m) = 0, Vé G [r,r + t] et Vm £ [0,g(l)]. 

Ainsi, N{t,m) est solution de pour te [t^, î^ + z] et m G [0,g(l)]. 

Montrons que N est unique. 
Soit Z{t,m) une autre solution de (fT8|) sur [t, r + t] x [0,.g(l)] : 

Z{t,m) = No{t,m) - J{Z){t,m), 

pour tout i G [t^, + et tout m G [0,(7(1)]. 
Alors, \/t G [T,r + t] et Vm G [0,3(1)], 

\N{t,m) - Z{t,m)\ = \ J{N){t,m) - J{Z){t,m)\, 

< al |7V(s,7r_(i_s)(m)) - Z(s,7r_(t_5)(m))|c;s. 

Pour i G [t, T + r] et m G [0,^(1)] nous posons : 

w{a) = |7V(cr,7r_(i_^)(m)) - Z(cr,7r_(t_^)(m))|, 



pour tout cr G [t, t]. 
Alors, 



w{a) ^ |J(Ar)(cr,7r_(i_^)(m)) - J(Z)(cr,7r_(t_„)(m))|, 
< j_ l^(s,7r_(<^_s) o7r_(t^^)(m)) 

-2'(s,7r_(^_^) o7r_(t_^)(TO))|dr. 
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Comme 



1"-(<T-s) ° '^-(t-a) = ^-{t-s) 



alors 



w{a) < al |iV(s,7r_(t_^)(m)) - Z(s,7r_(t_5)(TO))|dr, 



< al w{s)ds. 




D'après le lemme de Gronwall, Ver G [T,t], u)(cr) = 0, donc Ver G [t,t + t], w{a) — 0. Ainsi, 
VaG [T,r + r],Vme [0,5(1)], 



Donc N{t,m) = Z{t,m), G [r,T + r] et Vm G [0,5(1)]. 

Le problème (IT51) possède donc une unique solution sur [0,T + r] x [0,5(1)] qui est continue. 
Notons N* cette solution. 

En utilisant la suite {Nn)neN définie de la façon suivante : 
si t G [0,r + r] alors 7V„(t, m) = N*{t,m), Vm G [0,5(1)], et si t G [r + t,t + 2r] alors, 



nous montrons de la même manière que le problème (|18|) possède une unique solution, conti- 
nue, sur [0,T + 2r] x [6,5(1)]. Par la méthode des pas, l'existence et l'unicité d'une solution 
continue sur [0, -fcx)) x [0,5(1)] sont établies. □ 

Avant détablir le résultat d'unicité dépendant uniquement des cellules de petites matu- 
rités, nous montrons le résultat de positivité suivant. 

Proposition 5.1 Nous supposons que S{m) + V'{m) > 0, pour tout m G [0,5(1)]. Si, 
y{t,m) G [0,r]x [0,5(1)], N{t,m) > 0, alorsN{t,m) > 0, pourtoutt > t et tout m G [0,5(1)]. 

Preuve : Nous donnons une idée de la démonstration. Nous utilisons la méthode des ca- 
ractéristiques et la méthode des pas. 

Nous notons N = ip sur [0, t] x [0, 5(1)] et nous supposons (p > 0. Supposons alors qu'il existe 
(^0, nT-o) ë \t, t+t] x [0, 5(1)] tel que N{to, mo) < 0. Considérons alors l'équation différentielle 



A^(cr,7r_(t_^)(m)) = Z(cr, 7r_(t_^) (m)). 



Vm G [0,5(1)], 



{ 



No{t,m) = ip{T,n_(^t_^-^{m,))K{t -T,m) + G{N*){t,m), 
Nn{t,m) = iVo(i,m) - J(iV„_i)(i, m), 




et notons 



u[t) ^ N{t,m{t)), iG[r,T-Fr]. 



Alors 



u\t) 



{S{m{t)) + V'{m{t)) + (3{m{t), u{t)))u{t) 
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+2 j C{m{t), a)/3(A(a, m{t)),(p{t - a, A(a, m,{t))))ip{t - a, A(a, m{t)))da. 
Donc u'{to) > 0. En effet, u{tn) < et 

C(m(io), a)/3(A(a, m{to)), ip{tQ - a, A(a, TO(to))))</j(<o - a, A(a, m{to)))da > 0. 



Nous aboutissons à une contradiction. Donc N{t,m) > sur [r,T + r ] x [0,g(l)]. Par la 
méthode des pas, nous obtenons ce résultat pour t > r.D 

Nous allons à présent établir un résultat qui souligne l'influence des cellules de petites 
maturités (cellules souches) sur l'évolution de la population totale de cellules. 

Nous supposerons, dans la suite, qu'il existe C > tel que pour tout m € (0,g(l)] 

ds 

<C. (22) 



V{s 

Exemple 5.1 Si V{s) ^ as, avec a > 0, alors 



ds 



Vis] 



< C si et seulement si g'{0) > 0. 



Théorème 5.2 Nous supposons l'hypothèse vérifiée. 

Soit tq = sup„>Q ^ /^^ yfrj^ ■ Supposons que Ni{t, m) et N2{t, m) sont solutions de Iil8\) 
pour des données initiales ipi et ipi respectivement. Supposons que t > tq et qu 'il existe 
b G (0,/i-i(e-^)) tel que 

If lit, m) = if2{t,rn), pour m G [0,6] et t G [0, r] . 

Alors, il existe t > tel que 

Ni{t,m) = N2it,m), pour m e [0, g{l)] ett>t. 

Preuve : Remarquons que la condition t > tq implique, d'après le lemme lÏÏTTl que A(a, m) < 
m pour tout m G [0,^(1)] et tout a £ [Zi^]j 6t que h~^(e~-) < g(l). De plus, comme t > 0, 
h-\e-^) > 0. 

Montrons tout d'abord que Ni{t, m) — N2{t, m), pour m G [0, b] et t > 0. 
Soit m G [0, b]. Si t £ [0, r], le résultat est évident. 
Soit t >T. Rappelons que si N est solution du problème (fT5)) alors 

N{t,m) = (p{T,TT_çt--)im))K{t -T,m) + G{N){t,m) - J{N){t,m), 

pour tout t>Tei tout m G [0,.g(l)]. 
Comme i — t > 0, 

7r_(f_-)(TO) < m < b. 
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Donc 

Soit t £ [t, T + r] . Alors, s - a G [0, t - r] C [0, t] et 

A(a,7r_(t_^)(m)) < 7r_(t_^)(m) <m<b. 
Donc G{Ni){t,m) — G{N2){t,m). Par conséquent, 

liVi(t,m)-iV2(i,m)| = |J(iVi)(É,m)- J(iV2)(t,"i)| 

< a^y^ |Afi(s,7r_(i_s)("i)) - A^2(s,7r^(t-s)(TO))Ms- 

Par un raisonnement analogue à celui présenté dans la démonstration précédente et en utili- 
sant le lemme de Gronwall nous obtenons 

Ni{t,m) = N2{t,m), 

pour tout i G [t, r + r] et tout m G [0,6]. Par la méthode des pas, nous obtenons le résultat 
pour t>0. 

Soit a l'application définie de [0,5(1)] dans [0,5(1)] par 

a{m) = A(r, m). 

D'après le lemme lÏÏTTl a est une application croissante, atteignant son maximum en ^(l) et 
vérifiant aÇm) < m pour tout m G (0,5(1)]. a est donc inversible et nous notons A son 
inverse. Comme a(g{l)) = h^^{e^-), A est définie sur [0, h~^{e~-)]. Nous pouvons prolonger 
A par continuité sur [h~^{e~-), g{l)] en posant 

A(m) = g{l), pour tout m G [h~^{e~-), g{l)]. 

a étant strictement croissante et positive sur [0, 5(1)], A est strictement croissante sur [0, h~^{e^ 
Considérons la suite (6n)neN définie par 

bn+i = A(6„), pour n G N, 
bo = b. 

La suite (6„)„gN est croissante; en effet, Vn G N, 

bn = a{bn+i) = A(t, 6„+i) < 6„+i. 
Nous définissons aussi la suite {tn)neN par 

, , h{b^y 

(,„ = In ' 

c'est-à-dire 

I <n = 



hib) 

tn+i - t„ + r + In (^^gf^) , pour n G N, 
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La suite (t„)„eN est croissante. 
Montrons par récurrence que 

(iî„) Si m G [0, bn] et t > t„ alors Ni{t, m) = N2{t, m). 

Le résultat est vrai pour n = 0, nous venons de le montrer. 

Supposons {Hn) vérifiée. Soient m £ [0, 6^+1] eit> tn+i- Si t > tn+i alors t > tn + t. Nous 
pouvons réécrire les solutions Ni, i — 1,2, du problème (jlSp de la façon suivante : 



Ni{t,m) = Ni{tn + T,T:_(t-t„-r){m))K{t - tn ~T,m) - J{Ni){t,m) 
+ 2 ( ({n_(^t-s){m),a)N,{s ~ a,A{a,TT_(t-s){m)))x 

/3(A(a,7r_(f_5)(TO)),7Vi(s - a, A{a,T:_(^t-s){'m))))daj K{t ~ s,m)ds, 

pour t > t„ + T. 

Montrons que 7Vi(t„ + t, 7r_(t_t^_-) (m)) = A^2(in + î", 7i'_(t_t„_-) (m)). 
t > donc 

t — tn — T > In ' 

ou encore, 

g-(t-t„-r) < 



Kbn+l) ' 

D'où 

h{bn+i) h[bn+ij 

car h est croissante. Ainsi, 

En utilisant l'hypothèse de récurrence, nous obtenons le résultat souhaité. 
Enfin, comme s — a > tn et 

A(a, 7r_(t_5)(m)) < A(a, m) < A(t, m) < A(r, 6„+i) = 6„, 

nous obtenons 

\Ni{t,m) - N2it,m)\ = \J{N,){t,m) ~ J{N2){t,m)\. 

Par un raisonnement analogue au précédent, utilisant le lemme de Gronwall, nous vérifions 
l'hypothèse (iî„+i). 

Remarquons que, par définition de A, il existe M G N* tel que 

Bm < 5(1) = bM+i- 

Alors, pour tout m G [0, 6m+i] = [0,g(l)] et tout t > tM+i, 

Ni{t,m) = N2{t,m). 

Nous pouvons poser t = tM+i et le théorème est alors démontré. 
□ 

Nous déduisons de ce théorème le corollaire suivant : 
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Corollaire 5.1 Nous supposons l'hypothèse \2S]) vérifiée. 

Soit To = sup,„>o ( /m ''™^ ■ ^""^^ ^ ^ C°([0, t] X [0, 5(1)]) et N'^ la solution du problème 

associée à la donnée initiale ip. Supposons que t > Tq et qu'il existe b S {0, h~^{e~-)) 
tel que (p{t,m) — 0, pour m £ [0,5] et t [0,t]. Alors, il existe t > tel que N'-''{t,m) — 0, 
pour m e [0, 5(1)] et t >t. 



Ce résultat met en avant l'influence des cellules de petites maturités (cellules souches) sur 
la population future de cellules. Nous remarquons que dans le cas d'une absence de cellules 
souches dans le milieu initial (ce qui est le cas d'une anémie aplasique), l'extinction de la 
population après un temps fini t est inévitable, les cellules ne parviennent pas à survivre. 



6 Invariance — Application du théorème 15.21 

Dans toute cette partie, nous notons TV"^ la solution du problème (|18p pour une donnée 
initiale Lp G C"([0,r] x [0,51(1)]). Pour tout réel b e (0,g(l)], nous définissons l).]];, par 

VtieC°([0,r] X [0,5(1)]), \\u\\b :=sup{lM(t,m)l;(t,m) G [0,r] x [0,5]}. 

Nous supposons que 

/ inf (ô{m) + V'(m)) > 0, 

me[0,g(l)] 

et nous notons 

C:=sup{lC(m,a)l;(m,a) G [0,^(1)] x [G,t]}. 

Proposition 6.1 Supposons que Z(2(r — t)C + 1) < /. Alors, pour tout b G (G, g{l)] et toute 
fonction if G C°([0,r] x [0,^(1)]), 

Vi >T, VtoG [0,5], liV'^(i,m)l < 

Preuve : Soient b G (0,^(1)] et ip e C"([0,r] x [G,g(l)]). Pour i > G et m G [G, 6], nous 
définissons la suite {Nn)neN de la façon suivante : Vn G N, Nn{t,m) — (p{t,m) si t G [G,t], 
et, si t > r, 

NQ{t,m) = ip{T,TT_çt--){rn))K{t -T,m), 

Nn{t,m) = No{t,m) + G{Nn-^i){t,m) ~ J{Nn-i){t,m), n>l. 
Alors, Mt > T, Vm G [0,5], 

\Noit,m)\<e-'^'-~^\\ip\\b<\Mb. 

Montrons, par récurrence, que, pour tout t > r et tout m G [G, 6], 

VnGN, \Nnit,m)\ < \\ip\\b. 

Le résultat est vrai pour n ~ 0. Supposons qu'il est vrai au rang n, n > 0. Alors, pour tout 
t >T et tout m G [G, b], 

\Nn+i{t,m)\ < \Na{t,m)\ + \G{Nn)it,m)\ + \J{Nn)it,m)\. 
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Or, 

\G{N,,){t,m)\ < 2l{T~T)C\MtJ_ e-'^'-'Us, 

et ^ 

\J{Nn){t,m)\ < l\\ip\\bj_ e-"^*-'^ds. 

Donc ^ 

'(2(t-t)C + 1) 
En notant /i = j , nous obtenons 

|7V„+i(t,m)| < \\ip\\Je-'^'-^^ +tiil-e-'^'-^^) 



Comme /i < 1 par hypothèse, nous aboutissons à 

|A'"„+i(t,m)l < \\ip\\b- 

Nous obtenons le résultat attendu. 

D'après la démonstration du théorème 15.11 la suite (A^„)„gN converge vers N lorsque n 
tends vers +00. Par passage à la limite sur n nous obtenons finalement 

|iV^(t,m)| < \\ip\\b, yt>T,\Jm G [0,6]. 

La proposition est ainsi démontrée. □ 

Nous déduisons de ce résultat, le théorème suivant. 

Théorème 6.1 Si 1{2{t — t)Q + 1) < / alors, pour tout b G (0,5(1)] et toute fonction 
if G C°([0,t] X [0,g(l)]), il existe t > tel que 

yt>t, VTOG[0,g(l)], \N'P{t,m)\ <\\ip\\t. 

Preuve : Soient b G (0,5(1)] et 93 G C([0,t] x [6,5(1)]). Nous définissons, pour t G [0,t], 
l'application <^ par : 

_ J (p{t, m), si m G [0, 5], 
ip{t,b), si m G [b,gil)]. 



(p{t,m) 



D'après la proposition 16.11 pour tout t > r et tout m G [0,^(1)], \N'^{t,m)\ < ||<^||g(i). De 
plus, comme cp{t,m) = (p{t,m) sur [0,r] x [0,5], d'après le théorème 15.21 il existe ï > tel 
que N'^{t,m) — N'^{t,m) pour i > t et m G [0,^(1)]. Enfin, comme ||<^||g(i) — WfWb, nous 
obtenons 

Vt>l Vm G [0,5(1)], |iV'^(i,m)| < ll^llfc. 

□ 

Ce théorème met en avant le résultat suivant : si l'inégalité 1{2{t — t_)C, + 1) < I est vérifiée, 
alors la population de cellules souches domine la population totale, au moins à partir d'un 
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certain temps. Or l'inégalité 1{2{t — t)( + 1) < I est vérifiée, par exemple, si l est petit ou 
/ grand, c'est-à-dire, biologiquement, si peu de cellules sont introduites de la phase de repos 
vers la phase de prolifération (car l > f3{m, x) pour tous m G [0, 5(1)] et a; > 0) ou si le taux 
de mortalité dans la phase de repos est élevé. 
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